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In this study the ILC's capabilities for detecting WIMPs and measure their properties 
are investigated. The signal events are detected by associated production of Initial 
State Radiation (ISR). A model independent formulation of the signal cross section 
is used. The cross section is normahzed by inference from the observed abundance of 
cosmological Dark Matter (DM). The study is performed in full simulation of the ILDOO 
detector model. The prospects of determining the WIMP parameters individually and 
simultaneously are presented. 

1 Dark Matter and WIMPs 

New physics is required from both cosmol- 
ogy and the known shortcomings of the 
Standard Model (SM). The observed non- 
baryonic Cold Dark Matter (CDM) density 
in the universe can be accounted for by in- 
troduction of a new neutral particle, stable 
by virtue of a new conserved quantum num- 
ber. Given interactions of weak strength 
and masses in the range of 0.1-1 TeV these 
generic WIMPs naturally give the observed 
DM density as a thermal relic of the ex- 
pansion of the early universe. Several ex- 
tensions of the SM include such a WIMP 
candidate, e. g. SUSY with conserved R- 
Parity, where the LSP often is either the 
lightest neutralino Xi or the gravitino G. 
Direct WIMP pair production e+e^ — > XX 
should be possible at the International Lin- 
ear CoUidcr (ILC) [J for WIMP masses of 
less than half the CMS energy, M^ < y/s/2. 
However, since WIMP interactions arc only 
of weak strength, such events would be un- 
observable as long as no other detectable 
particle is associated with the pair produc- 
tion process. ISR provides such a detectable particle. A single photon radiated off the initial 
state delivers a clean signal, while the properties of the photon spectrum give information 
on the DM candidate. This detection channel does not only give additional information on 




Figure 1: Quadrant of the ILD detector. 
Starting from the origin of the coordinate sys- 
tem the Time Projection Chamber (yellow) 
and the calorimetric system of ECAL (green) 
and HCAL (ocher) are shown. 
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Figure 2: Schematic derivation of the e+e^ — > XXI production cross section. Except for 
some generic assumptions on the WIMP nature no presumptions on the interactions involved 
are imposed, indicated with the question marks. 



the new particles, but it might as well be the only kinematically allowed process at the ILC. 
This is the case when only one new particle exists in the energy range of the ILC. In super- 
symmetric scenarios the other SUSY masses might be just above the production threshold, 
preventing the production of the LSP via the decay of heavier particles. Without the need 
for heavier particles to be produced, ISR searches extend the ILC reach on WIMP masses to 
M^ < 250 GeV, just below half the CMS energy. A model independent description of ISR 
associated WIMP pair production can be derived as in [2] . Under the assumption of only 
one new particle responsible for the CDM the annihilation cross section of a WIMP pair into 
an e~^e~ pair can be estimated from the observed DM relic density (Fig.[2][a)). Exploitation 
of crossing symmetry yields the cross section for the inverse process (Fig. [DJb)). For soft or 
collinear photons, factorization theorems can be applied to get the observable signal process 
depicted in Figure Wl,c) . It can be shown that the soft/collinear approximation holds also 
for photons with larger polar angles [2]. The resulting differential production cross section 
as a function of the emitted photon energy is given by: 



da 
ax 



Pp)2'''"{2S^ + ir 1 



(l-x)s 



l/2+,/o 



(1) 



where x = 2Ey/^/s and Ke{Pe,Pp) is the polarization dependent squared coupling to the 
electrons in the initial state. The parameter Jq is the angular momentum of the domi- 
nant partial wave in the production process. Depending on the nature of the interaction 
Jo is related to the spin of the exchange particle. Since all known models predict either 
s-wave (Jo = 0) and p-wave (Jo = 1) WIMP annihilation only these cases are considered 
here. Figure |3] show the cross section as function of the recoil mass M^^^^^^ = s — 2^/sE^ 
for (a) WIMP masses of 120, 160 and 200 GcV and (b) for a 200 GeV WIMP with Jo = 
and Jo = 1. The cut-off in the spectrum is related to the WIMP mass and the shape at the 
threshold provides information on the partial wave quantum number. In order to measure 
the energy cut-off and the shape, a very good energy resolution for photon detection is de- 
manded from the detector and the reconstruction algorithms. 

This analysis is performed with a full simulation of the International Large Detector con- 
cept (ILD) [3]. A quadrant of the ILD is shown in Figure [TJ The ILC's capability of 
delivering polarized beams can be used to significantly suppress SM backgrounds, and in- 
crease possible signals depending of the structure of the involved interactions. The role of 
beam polarization at the ILC is discussed in |3]. 
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3: Differential xxi production cross section in terms of the recoil mass for (a) different 
masses and (b) Jq = and Jq = 1 (in arbitrary units). 



X 



2 ILD detector simulation and plioton reconstruction 



Single plioton events provide the means for 
a fast validation of the detector model's im- 
plementation in the simulation framework. 
Some key parameters like the overall mate- 
rial budget and the intrinsic energy resolu- 
tion of the calorimetric system for electro- 
magnetic processes can be determined from 
the analysis of a large simulated sample of 
photons shot randomly into the detector. 
Secondly some aspects of the employed re- 
construction algorithms can be studied us- 
ing single photon events, for example the 
quality of pattern recognition for electro- 
magnetic clusters, addressed in Subsec. 12.21 
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2.1 Material budget 



Figure 4: The material budget of the Mokka 
implementation of the ILDOO detector model 
tracking system in terms of the radiation 
length Xq as a function of the polar angle 8. 
For comparison the budget as described in the 
Lol [3] is plotted in the blue histogram. 



A good understanding of the material bud- 
get is important for the evaluation of the 
systematic effects from photon conversions 
on the event reconstruction. The angular 
dependence of the material budget that par- 
ticles originating from the interaction point 
are exposed to on their way out through the 
tracking subsystems (vertex detector, TPC, 
intermediate trackers and support struc- 
tures) can be evaluated in the detector simulation by counting the number A^gen of photons 
generated and the number Nunconv. of photons not converted before reaching the ECAL 
face. The required information is given in the output of the Mokka implementation of the 
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ILDOO. The material budget in terms of the radiation length is given by: 

Figure|4]shows the amount of material determ.ined for the ILDOO detector model implemented 
in MoKKA v-06-07-p01 The results are based on 10^ single photon events from the SM 1/1/7 
process. For reference the material budget for the tracking system as quoted in the ILD 
Lol [3] is shown in the underlying histogram. The two most visible features are the increased 
amount material at cos 9 ~ 0.0 which is due to the dividing central cathode in the TPC, and 
the additional material cos0 > 0.8 in the forward region is related to the TPC endplate. 

2.2 Photon reconstruction 

The PFlow paradigm requires that every particle partaking in an event is to be recon- 
structed individually. For single photon searches especially the performance of the recon- 
struction algorithms in finding electromagnetic clusters in the ECAL is required to be very 
good. Figure [5] depicts the average ratio Nrec/Ngen of reconstructed photon candidates per 
generated photon as a function of the generated photon energy (a) and generated photon 
polar angle cos (8) (b). For event reconstruction the Marlin framework vOO-10-04 and the 
PandoraPFlow [5] algorithm v-03-01 have been used. A photon is considered as recon- 
structed when a PFlow photon object is found within 0.1 rad (seen from the IP) from a 
generated one. Additionally the photon object is required to be related to the generated 
photon on MC level. The ratio continuously rises from one at low photon energies to about 
three at E^ = 250 GeV. The angular dependence displays a steep rise at polar angles of 
cos ~ 0.8. This angular region corresponds to the insensitive region between the ECAL 
barrel and ECAL cndcap where pattern recognition is difficult due to the complex geom- 
etry (see Fig. [1]). For the analysis presented in the next section this effect is countered 
by combining the "split" photon clusters into higher level photon objects. Photon clusters 
that are less apart than 0.04 rad are considered as being fractured remnants from a single 
incident photon. The result of the merging procedure is shown in Figure [HI where the ratio 
deviates from unity on the sub-percentage level. 

3 Determination of WIMP parameters 

WIMP production with associated ISR leads to an excess in the photon spectrum at energies 
below 

El - Ml 
E, = — ^^- (3) 

The cross sections are in the order of a few fb, depending on the values of the model 
parameters Ke{Pe,Pp), Jo, M^ and the WIMP spin S^, resulting in an S/B ratio of 10~^ — 
— 10^^ to the irreducible lyiyy background. The parameter Ke{Pe,Pp) can be deconstructed 
according to the four possible e^e~ helicity configurations: 

KeiPe,Pp) = J(l + Pe)[(l + PpHeRPL) + (1 - Pe+)<eHPE)] 

+ i(l - Pe)[(l + Pp)>i{eLPL) + (1 - Pe+HeLPR)] (4) 
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Figure 5: Performance of the photon reconstruction in terms of the ratio Nrec/Ngen of 
reconstructed per generated photons as function of photon energy (a) and polar angle (b) . 



The WIMPs might couple only to right-handed electrons and left-handed positrons (^(6^,^^), 
helicity conserving). The couplings might in addition also conserve parity (k^ch^pl) = 
i^{^l,Pr))- In both cases the S/B ratio can be enhanced using polarized beams. The results 
presented here are based on the full SM i'iyy{N)^ {N = 0,1,2) background simulated for 
the ILDOO detector model implemented in MOKKA v-06-07-p01. For event generation the 
Whizard event generator [6] is used. Event reconstruction is done within the Marlin 
framework vOO-10-04 using the PandoraPFlow algorithm v-03-01. Each event contains 
up to two additional ISR photons from the electron and positron. In total an integrated 
luminosity of £ = 250 fb~^ is generated and simulated. The signal itself is not generated and 
simulated, but is obtained by assignment of weights to the dominant vv^ background, where 
the weights are given by the ratio of the cross sections for neutrino and WIMP pair produc- 
tion, w = (^xxil^vv-y evaluated for E^ and 6 of the detected ISR photon. For the signal to 
be statistically independent of the background only half of the simulated SM background 
is used for the weighted signal, while the other half serves as signal free background. An 
additional weight is assigned to each event to adjust for the demanded luminosity. The event 
selection requires a single photon in the detector and no activity in the tracking system. 

3.1 2ct reach on k, 

The reach on 2ct level of the ILC on the parameter k as function of the WIMP mass is shown 
in Figure [71[a) . The plot corresponds to an integrated luminosity of 200 fb~^. The WIMP 
parameters are chosen as 5 = 0, Jo = 1 and the WIMP couplings to electrons conserve 
helicity and parity. The contour lines give the reach on the parameter k as a function of the 
WIMP mass. The area above the curves is accessible. For each WIMP mass a signal region 
is defined in the energy distribution of the emitted photons, Emin < E.^ < Emax- The lower 
cut on the photon energies ensures that the recoiling WIMP system is non-relativistic. This 
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Figure 6: Ratio of Nrec/Ngen for photons after merging of fractured photon objects. 



has to be enforced since only s-wave and p-wave production is assumed. The upper cut is 
given by the kinematic limit (Eq. [3]). This cut serves to increase the signal to background 
ratio. The result is plotted for different beam polarizations: 

• (Pe/Pp) = (0.0/0.0) (red, solid) 

• (Pe/Pp) = (0.8/0.0) (blue, wide dashed) 

• (Pe/Pp) == (0.8/ - 0.3) (green, dashed) 

• (Pe/Pp) = (0.8/ - 0.6) (black, dotted) 

In this scenario and the given luminosity, no value of k is accessible for unpolarized beams, 
hence no WIMP signal can be distinguished from the SM background on 2a level (the 
corresponding contour line coincides with k = 1). With increasing beam polarization the 
reach is extended for WIMP masses between 160 and 200 GeV down to values of 0.6. The 
loss of sensitivity to higher masses is due to the kinematic limit, while for lower masses the 
signal is increasingly reduced by the lower cut on the photon energies. Figure [ZKb) shows 
the reach for the case where the WIMPs (S — 1, Jo = 1), couple only to right handed 
electrons and left handed positrons. Even for unpolarized beams the reach covers large 
areas of the M — k plane for masses up to 200 GeV. With increasing beam polarizations 
the SM background is suppressed and the signal enhanced, resulting in a large increase of 
accessible values down to k < 0.1. 

3.2 Mass measurement 

This section on the WIMP mass measurement refers to results from 2007 [7] . The simulated 
detector model is LDCPrime_02Sc. The mass of a WIMP candidate can be determined by 
comparing the signal spectrum to template spectra of different WIMP mass hypotheses 
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Figure 7: 2a reach for 200 fb~^ on Ke{Pe,Pp) as function of the WIMP mass. The results 
are given for (a) spin WIMPs with hehcity and parity conserving couplings and (b) spin 
1 WIMPs with couplings conserving only helicity. The contour lines are given for several 
beam polarizations, see text. 

but with otherwise identical parameters. Figure [Sfa) shows the x^ foi' the template fits as 
a function of the template mass. The signal is a 180 GeV spin 1 WIMP with couplings 
conserving helicity and parity. The other parameters are set to k = 0.3 and Jq = 1. 
The simulated data is equivalent to an integrated luminosity of 500 fb~^. Three different 
polarization configurations are considered: 

• (Pe/Pp) = (0.0/0.0) (solid) 

• (Pe/Pp) = (0.8/0.0) (dotted) 

• (Pe/Pp) = (0.8/ - 0.6) (dashed) 

Depending on the beam polarization the mass can be determined to: 

• M ^ 181.5 ± 3.0 GeV ((Pe/Pp) == (0.0/0.0)) 

• M = 180.5 ± 0.5 GeV ((Pe/Pp) = (0.8/ - 0.6)) 

In this case due to the preferable choice of the WIMP couplings highly polarized beams 
increase the mass resolution by a factor of six. 

3.3 Simultaneous determination of mass and Jq 

As pointed out in Section [T] the shape of the photon distribution at the threshold provides 
information on the dominant partial wave of the production process. The energy distribution 
of ISR photons from the SM background and the WIMP production can be fitted with the 
WIMP cross section of formula [TJ leaving the WIMP mass and the partial wave Jq as free 
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Figure 8: x^ of comparison of 180 GcV spin 1 WIMP {k ~ 0.3, Jo = 1) signal distribution 
with signal mass templates for three different polarization configurations, as explained in 
the text. 



parameters in the fit. For C = 100 fb ^ {{Pe/Pp) = (0.8/ — 0.6)) this method recovers 
a WIMP mass of Mx = 173.8 ± 5.1 GeV and Jo = 1.26 ± 0.34 for input parameters of 
Mx = 180 GeV and Jo = 1 respectively. The remaining input parameters are k = 1.0, S = 1 
and the couplings conserve helicity only. Although the luminosity is low compared with the 
ILC design, the s-wave production can be excluded on 3ct level, due to the assumed WIMP 
couplings and beam polarization. The simultaneous mass fit performs significantly worse 
than the dedicated template fit of the previous section. Since the mass can be determined 
independently from Jo, a constraint could be imposed on it in the simultaneous fit. The 
expected improvement on the Jo measurement has to be quantified. 



4 Conclusions 



The study presented here shows that there is a good chance of a model independent detection 
of WIMPs at the ILC. If the WIMP couplings to electrons are not too small, their masses can 
be measured on the percentage level. Both the detection sensitivity and the mass resolution 
is enhanced when polarized beams are assumed. First results indicate that even the more 
involved measurement of the partial wave quantum number seems to be possible, given a 
good understanding of the SM background. 
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